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Abstract 

We study the one-loop renormalization and evolution of the couplings in scalar field 
theories of the Lifshitz type, i.e. with different scaling in space and time. These theories 
are unitary and renormalizable, thanks to higher spatial derivative terms that modify 
the particle propagator at high energies, but at the expense of explicitly breaking Lorentz 
symmetry. We study if and under what conditions the Lorentz symmetry can be considered 
as emergent at low energies by studying the RG evolution of the "speed of light" coupling 

and, for more than one field, of 5c 2 = — c^ 2 in simple models. We find that in the UV 
both c 2 and 5c 2 generally flow logarithmically with the energy scale. A logarithmic running 
of c 2 persists also at low-energies, if 5c 2 ^ in the UV. As a result, Lorentz symmetry is 
not recovered at low energies with the accuracy needed to withstand basic experimental 
constraints, unless all the Lorentz breaking terms, including 5c 2 , are unnaturally fine-tuned 
to extremely small values in the UV. We expect that the considerations of this paper will 
apply to any generic theory of Lifshitz type, including a recently proposed quantum theory 
of gravity by Hofava. 



1 Introduction 

Recently, there has been an increasing interest in non-relativistic quantum field theories 
where Lorentz invar iance is explicitly broken at high energies and hopefully recovered at 
low energies. In particular, in [U O [3] (see also [4j), general gauge theories, including 
non-relativistic extensions of the Standard Model, were proposed and investigated, while 
in [U [6] similar constructions were implemented in Yang-Mills theories in 4+1 space-time 
dimensions and membrane theory. The same type of construction was then extended to 
four-dimensional quantum gravity in [7J where it was suggested that the resulting theory 
may provide a candidate for a renormalizable and unitary quantum theory of gravity which 
flows in the infrared (IR) to Einstein theory. 

The ultraviolet (UV) behavior of all these theories is substantially ameliorated by the 
presence of higher derivative (in the spatial directions only) quadratic terms that improve 
the UV behavior of the particle propagator, without introducing ghost-like degrees of 
freedom that in Lorentz invariant higher derivative theories typically spoil the unitarity 
of the theory. The proposed theories are of Lifshitz type [TO]- In the UV, they exhibit, 
at the classical level, an anisotropic scaling symmetry under which time and space scale 
differently: t = X z t', x = Xx' , where z is the critical exponent, equal to one in a Lorentz 
invariant theory. The renormalizability properties of these theories have been extensively 
studied in [TJ[2] for scalar, fermion and gauge theories. The usual power-counting argument 
for the renormalizability of a theory does not strictly hold anymore, but it is essentially 
still valid, provided one substitutes the standard scaling dimensions of the operators by 
their "weighted scaling dimensions" [1], i.e. by the dimensions implied by the assignment 
[x]w = —1 and [t] w = —z. Lifshitz-type theories exhibit at least two qualitative different 
energy regimes, set by the scale (denoted by A) of the higher derivative operators. We will 
generally denote as UV regime the energy range E ^> A, where the theory is manifestly 
non-Lorentz invariant. This is the proper "Lifshitz" regime, where the effective scaling 
dimensions of the operators are the weighted ones. We denote as IR the range E <C A, 
where the theory is expected to smoothly reach the "standard" regime, where the operators 
are classified by their standard scaling dimensions. Weighted relevant operators break the 
anisotropic scaling symmetry explicitly and, at low energies, the theory is expected to flow 
to the Lorentz invariant theory with z = 1. This is, however, a non-trivial (and obviously 
crucial) point, since there is no dynamical principle for which Lorentz symmetry should 
emerge in the IR. As a matter of fact, we find that Lorentz invariance is recovered in 
the IR only if an unnatural fine tuning of the parameters of the theory ensure that all 
sources of Lorentz violation in the theory are tiny and below the rather stringent, presently 
known experimental bounds. These issues, expected on general grounds [11], are explored 
by performing a one-loop calculation in concrete, simple, models, which also clarify the 

1 See also 8,9 for a (partial) list of further works that investigate the proposal of [7]- 



2 



IR-UV structure of Lifshitz-type theories. 

More specifically, the purpose of this paper is to study the Renormalization Group (RG) 
evolution at one-loop level of simple scalar field theories of Lifshitz type. In particular, we 
will calculate the one-loop beta-functions of the weighted marginal operators in the theory 
and solve the corresponding equations to study the evolution of the associated couplings 
gi. Once this is performed, we will focus our attention on a particular weighted relevant 
operator, c^(d(j)) 2 , and study the RG evolution of the "speed of light" parameter for <ft 
(physically, in the low energy theory represents the maximal speed for cf) particles). In 
order to keep the technical analysis as simple as possible, we will mostly consider scalar 
field theories in higher dimensions with z = 2, namely a <^> 6 theory in D = 4 spatial 
dimensions, with quartic derivative couplings as well, and a (ft 3 theory in D = 10 spatial 
dimensions. These two theories are among the simplest theories which are i) of Lifshitz 
type, ii) their /3-functions are respectively positive and negative in the UV, iii) c 2 , has a 
non-trivial running already at one-loop level. These theories are obviously toy laboratories 
(in particular, the </> 3 theory is not even stable), yet they manifest, in a simple context, the 
main features that more complicated and "realistic" models of this sort should exhibit. In 
both theories, typically shows a logarithmic running in the range E 3> A, where A is 
the high-energy scale where the theory is anisotropic: 



c 2 (E) = c 2 (E ) 



l + /log(-|- 



(1.1) 



here an 0(1) particle-dependent constant and E A a reference scale in the UV 
range. In eq. (|l.ip . we schematically denote by / the radiative coefficient governing the 
RG flow, which depends on the coupling constants of the weighted marginal operators. It 
is also important to investigate what happens in the presence of more than one field, and 
particularly if and under what conditions their "speed of light" parameters converge to 
the same value. To address this question, we have also studied the RG evolution of the 
difference 5c 2 = c 2 ^ — c 2 2 . Under the assumption that Sc 2 <C 1, one schematically finds, 
for E > A, 



Sc 2 (E) = 5c 2 (Eq) 



"1 n s 



+ Sg(E ) 



l + /"log(^ 



U, (1.2) 



where 5g are small perturbations around some fixed-point solutions of the RG evolution. 
Eqs. (|l,ip and (jl.2p summarize the quantum evolution of c 2 and 5c 2 in the UV regime, as 
given by the marginal couplings. 

After having studied the UV, we move on to analyze the IR regime E <C A. We will 
see that A is the characteristic scale below which Lifshitz theories turn into "standard 
theories". More precisely, we will explicitly show that the RG evolution of all weighted 
marginal couplings is essentially frozen below A, in complete analogy to the decoupling of 
a massive particle in a standard quantum field theory. The key point is, of course, that in 
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the IR the relevant propagator term is the usual one, quadratic in the momentum, while 
the higher derivative terms can be neglected (it is however a delicate point, given that the 
theory with the usual quadratic propagator is non-renormalizable) . Taking into account 
only the effect of weighted marginal couplings, for E <C A we find 



which shows that, for sufficiently high A, the IR effect of the weighted marginal couplings 
can be neglected. This is expected, since in the IR the usual classification of operators in 
terms of canonical rather than weighted dimensions holds. What is marginal in the UV 
becomes then irrelevant in the IR. We will explicitly show how the /^-functions smoothly 
change their behavior going from the UV to the IR by computing them in a momentum 
subtraction renormalization scheme, where all the decoupling effects are manifest. 

However, care has now to be paid for the weighted relevant operators which become 
standard marginal in the IR, since they can efficiently mediate the UV Lorentz violation 
to the low-energy theory. Indeed, we will show, by explicitly working out a toy example 
in 3+1 space-time dimensions, that a logarithmic running like eq. ljl.ip (with Eq replaced 
now by A) still holds in the IR, with / depending now on the (standard) relevant couplings 
and being proportional to any Lorentz symmetry breaking coefficient of the low-energy 
effective theory, remnant of the Lifshitz-like nature of the UV completion. 

In general, then, Lorentz symmetry is not emergent in the IR in theories of Lifshitz 
type. Recovering Lorentz symmetry would require some dynamical principle keeping all 
sources of Lorentz violation sufficiently small. The experimental bounds on 5c 2 for ordinary 
particles are of the order of l(H 21 - 23 ) [12], which give an idea of the order of magnitude 
of the fine-tuning that is needed. An indirect bound on 5c 2 for any charged particle is 
implicitly given by eq. ljl.ip . In the case of photons, for instance, an experimental constraint 
on the energy dependence of c 2 by the FERMI experiment [13] gives, taking n 7 ~ 1 in 
eq. ljl.ip . the following rough constraint on /: 



Modulo a loop-factor and coefficients of order one, the bound (jl.4p can be seen as a 
bound on Sc 2 for any charged particle. We expect that these fine-tuning problems will 
affect all generic quantum field theories of Lifshitz type, in particular, the non-relativistic 
standard model of [3] and the proposed quantum gravity theory of [7]. In the latter case, 
after coupling the theory to matter, the problems mentioned above will reappear for the 
Standard Model particles, where parameters like 5c 2 have tight experimental constraints. 

The plan of this paper is as follows. In section 2 we briefly review the main properties 
of the Lifshitz-like theories, using, for the sake of illustration, a scalar field theory in 3+1 
dimensions. In section 3.1 we study the one- loop renormalization of a single scalar field 




(1.3) 



l/l < 10 



16 



(1.4) 
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</> 6 theory with derivative cf) 4 interactions in D = 4 spatial dimensions; in section 3.2 this 
analysis is extended to the case of two coupled fields. In section 4.1 we study the one-loop 
renormalization of a single scalar field c/) 3 in D = 10; the two-field case is dealt with in 
section 4.2. The analysis in sections 3 and 4 are performed using the minimal subtraction 
(MS) scheme, suitable for finding the /^-functions in the UV regime (E 3> A). In section 
5 we study the IR behavior (E <C A) of Lifshitz-like theories by using the momentum 
subtraction scheme. After reviewing in section 5.1 the IR behavior of the (3 function in 
conventional 4> A - theory, in section 5.2 we show that the /3-functions produced by weighted 
marginal couplings go to zero for E <C A. In section 5.3 we discuss the contribution of 
the weighted relevant operators (marginal in the standard sense) on Lorentz symmetry 
breaking effects in the IR. We will see that the presence of a non-zero 5c, inherited from 
the UV, induces a running in c 2 at low-energies. The effect is general and we expect that 
it should apply to any low-energy field theory (i.e. not only to theories of Lifshitz type) 
perturbed by Lorentz symmetry breaking terms. In section 6 we conclude. 

2 Renormalizable Lifshitz-like scalar field theories 

Unconventional scalar field theories of the Lifshitz type, with higher derivative interactions 
and higher derivative quadratic terms, have been extensively studied in [TJ. Here we briefly 
review some aspects of the construction that will be useful in what follows and refer the 
reader to [1] for more details. As mentioned in the introduction, the key point of the 
whole construction is to break Lorentz invariance, so that one is allowed to introduce 
higher derivative terms in the spatial derivatives and quadratic in the fields, without 
necessarily introducing the dangerous higher time derivative terms that would lead to 
violations of unitarity. In doing so, the UV behavior of the propagator is improved and 
theories otherwise non-renormalizable become effectively renormalizable. A useful guiding 
principle to easily classify and identify the renormalizable theories in this enlarged set-up 
is achieved by demanding an invariance under "anisotropic"!^] scale transformations: 

t = X z t', x^Xx 1 ', (l){x\t) = X E ^ 1 (j)'(x i ',t'), (2.1) 

where i = 1,...,D parametrizes the spatial directions. The parameter z is known as 
"critical exponent" and, when it equals one, the transformations (|2.1|) reduce to the usual 
Lorentz invariant scale transformations. According to eq. (|2.1j) . we can assign to the coor- 
dinates and to the fields a "weighted" scaling dimension as follows: 

[t] w = -z, [x i ] w = -l, [4>] w = ^l. (2.2) 

2 The word "anisotropic" arises from condensed-matter physics. In all the instances we consider, we 
assume to be in a so-called "preferred frame" |12j where spatial 50(3) rotations, translations and time- 
reversal are unbroken symmetries. 
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It is straightforward to see that at the quadratic level, modulo total derivative terms, the 
Lagrangian for a single scalar field, invariant under (|2.1I) . reads 

C quad = \ft - ^(Li) (^) 2 , (2-3) 

with a 2 being a dimensionless coupling and A a high-energy scale parametrizing the 
strength of the higher derivative operator. Due to the improved UV behavior of the 
propagator resulting from (|2.3p when z > 1, the usual power-counting argument for the 
renormalizability of a theory is no longer applicable. The required modification is obtained 
by substituting the scaling dimensions of the operators by their "weighted scaling dimen- 
sions" In other words, a theory is renormalizable if all the operators Oi appearing 
in the Lagrangian have weighted scaling dimensions [Oi] w (not to be confused with the 
standard scaling dimensions [Oi]) which are not greater than z + D. Thus, the second 
term in (12.3H . although manifestly irrelevant in the standard sense, behaves (and should 
be considered) as a marginal operator in this theory. 

It is useful to illustrate this construction with a specific simple example, namely a 
scalar field in 3 + 1 space-time dimensions (D = 3) and z = 2. For simplicity, we also 
impose a Z 2 discrete symmetry (j) —* —(p. The most general renormalizable Lagrangian, 
invariant under the transformations (|2.ip . is given by 

Cr = Ifi - ^-(A0) 2 - -h-id^fcf - -^<f> 10 , (2.4) 
2 Y 2A 2V Y> 48A 4V lYJ Y 10!A 6 ^ ' v ; 

where A = didi is the Laplace operator in the spatial directions. All the operators ap- 
pearing in (|2.4p are weighted marginal. The renormalizability properties of the theory are 
not changed if the scaling symmetry (|2.1|) is softly broken by adding weighted relevant 
operators. They are given by 

c = - ^mt + e (2 , i+2 >„-^ 2 " +2 + ■ («) 

n=l v ' 

so that the final Lagrangian is the sum of eqs. (|2.4p and (|2.5p . In conventional scalar 
field theories in 3 + 1 dimensions, the interactions appearing in (|2.4p would be non- 
renormalizable. What renders the theory renormalizable - and the interactions in (|2.4p 
weighted marginal - is the modification of the propagator, which in momentum space now 
reads 

^fcg-c 2 * 2 --^* 4 -™ 2 )" 1 , (2.6) 

with k 2 = kiki. The 1/k 4 high -energy behavior of the propagator leads to an improvement 
of the ultraviolet behavior of the theory. As a result, if no coupling of higher dimension is 

3 Sometimes in the literature the weighted scaling dimension is introduced as the standard scaling 
dimension in some non-standard natural units. Although there is nothing wrong in doing so, we prefer to 
distinguish between [Oi] w from [Oi] and use the usual natural units. 
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added, the theory is power-counting renormalizable. In the UV the RG evolution of the 
weighted relevant parameters, such as c 2 or g\ in eq. (|2.5p . will be governed by the RG 
evolution of a combination of the weighted marginal couplings /12, <?4 and a 2 . As we will 
see, in the IR the Lifshitz-type theory turns into a low-energy effective theory, where the 
weighted marginal couplings turn back into standard irrelevant ones and do not effectively 
run anymore. In this regime, if no Lorentz violating parameter appears in the Lagrangian 
terms with standard dimension < 4, then we effectively recover Lorentz symmetry in the 
IR, which protects c 2 from any possible running. On the other hand, if some Lorentz 
violating parameter is left (like e.g. 5c = — C0 2 7^ in a two-field model), it will still 
generically induce a running of c 2 , governed now by standard marginal couplings. 

For the 3+1 dimensional model defined by eqs. (|2.4p and (|2.5p . there is no renormal- 
ization of the couplings at one-loop level, due to the fact that the vertices in (|2.4p involve 
at least six fields. For this reason, in what follows we will consider higher dimensional 
scalar field theories, for which the weighted marginal vertices contain less fields and, as we 
will see, there is a non-trivial renormalization of the couplings already at one-loop level. 

3 UV behavior: a model with z — 2, D = 4 
3.1 One scalar particle 

We look for a weighted renormalizable scalar field theory with a non-trivial renormalization 
of the (di(p) 2 operator at one-loop level. A simple quantum field theory of this sort is 
obtained in 4 + 1 space-time dimensions with anisotropic scaling z = 2. The most general 
renormalizable Lagrangian, up to total derivative terms and including all possible weighted 
marginal and relevant operators, is 

C = Ij? - 4^(A0) 2 - ~(^) 2 - —<f> 2 - 4t^ 4 - 4^(W) 2 - 4n<f ■ (3-1) 
2 2A 2V ; 2 y ' 2 4!A 4A 3 K 1 6!A 4 

All couplings appearing in C, but the mass m and A, are dimensionless. In order to 
reduce the number of operators, we have imposed on £ a discrete Z2 symmetry under 
which <j> — > —(p. For simplicity, in the following we set A = 1. 

Our first aim is to renormalize the theory at one-loop level and study the RG flows of 
the weighted marginal couplings g and k in the UV. The coupling a 2 , although weighted 
marginal, is one-loop finite, since there is no way to extract four powers of external spa- 
tial momentum from the tadpole graph given by the quartic couplings appearing in C. 
Similarly, the wave function renormalization of <j> is trivial at one- loop level, Z = 1 + 0(2- 
loops). Once the RG flows for g and k are solved, we will study the evolution of the 
weighted relevant coupling c 2 @ We regularize the theory using a variant of dimensional 

4 Strictly speaking, the RG evolution of c 2 , as determined in a physical scheme, starts at two-loop order, 
even in presence of the quartic derivative interaction, since the momentum-dependence of the one-loop 
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regularization applied only to the spatial directions (D = 4— e) and thus renormalize using 
a minimal subtraction scheme where only the poles in 1/e are subtracted, with no finite 
term. 




(6) 




Figure 1: One- loop graphs contributing to the renormalization of the 4> 6 vertex. All external 
momenta are vanishing. 



The superficial degree of divergence of a graph is easily computed by looking at the 
weighted scaling of a graph. The one-loop corrections to the coupling k come from a one- 
loop graph with 3 insertions of the coupling g and from another graph with one insertion 
of k and g. See figure [TJ Due to some unusual properties of these theories, we report, in 
detail, the computation of the divergence term of the graph (a) in fig. [TJ A divergence 
can only arise when all the momenta of the vertex are taken in the internal lines, so that 
we can set to zero all external momenta pf 



£\3 



(a) = {-ig^ € ) 



15 f dtfdPq 6 3 

~2j WY^ q G{q ' q 



(3.2) 



where 15/2 is a geometrical factor taking into account all possible channels and 



]q — a 2 q 4 — c 2 q 2 — m 2 



(3.3) 



is the propagator for 0. Here and in the following q 2 = qiqi. After Wick rotating q° 
(q° = iq 5 ), we can rewrite (a) as 



(a) 



15^V e d 2 



4 

I5ig 3 



(dm- r j 
I3 i + finite 



(27T) 



D+l 



(3.4) 



graph (which is a tadpole) is trivial. However, we can still define a running c 2 coupling by adding a 
fictitious momentum in the loop, seen as the momentum carried by the composite operator (di4>) 2 . This is 
a standard trick. See e.g. [14] for the completely analogous case of the one-loop RG evolution of the mass 
parameter in the usual cf> 4 theory. 



8 



where 



f°°, j+ i [ dq 5 d D g 

I daa JWY^ q 



in e -a(ql+a 2 q A +m 2 ) 



(3.5) 



In writing the last equality of eq. (|3.4p . we have expanded the ac 2 q 2 term in the exponen- 
tial, since insertions of these terms lower the divergence of the graph. It is straightforward 
to check that the divergence arises only from the leading, c 2 -independent term. The inte- 



gral (|3,5p is easily done by going to radial coordinates and changing variables aa 2 q A 
Performing the integrals we get 



'"J 



4(2vr) 



D+l 



[m 



(a 2 Y 



p+2n (D + 2n\ /6 + 4j - D - 2n 



(3.6) 



where tt D = 2tt d / 2 /T(D/2) is the area of the S D sphere. Using the same techniques, we 
can compute the graph (b) as well. By denoting —5^^/61 the Lagrangian counterterm 
canceling the divergences coming from the graphs (a) + (b), we then get 



45# 3 / 4 15#/cZ 4 



32a 5 



+ 



8a 3 



(3.7) 



where we find convenient to express the result in terms of the usual loop factor for Lorentz 
invariant theories in D space-time dimensions, Id = £~lf)/(2ir) D . Similar manipulations 
allow to compute 5 g , the coefficient of the counterterm — 5 g ((j)di(j)) 2 /A: 



_ 3g 2 k 1 
9 ~ 8a 3 e - 

From eqs. (|3.7|) and (|3.8|) we obtain the one-loop (3 functions for g and k: 

15Z 4 gk 45/ 4 g 3 



(3.8) 



k = f3 k 
9 = 



8 a 3 
a 3 



32 a 5 



@9 - q „3 



(3.9) 



where a dot stands for a derivative with respect to t = ln/i//io and ^0 3* 1 is a given UV 
reference scale. Note that the effective couplings of the theory are 



9 

9 = -3 



k 



(3.10) 



The solutions of the RG equations (|3.9|) are 

9o 



9(t) 
k{t) 



1 

ko 



3Ugo 4. 
8 L 



9o 



+ 



5g 2 (git) 



9o 



9o 



(3.11) 
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with <?o = <?(0), fco = k(0). Since a Landau pole appears at the scale 

E p oie = We 5 ® , (3.12) 

the range in which eqs, (|3.1ip are reliable is 1 <C E <C E po i e . 

Having found the RG evolution of the weighted marginal couplings k and g, we can 

now go on and study the evolution of the weighted relevant coupling c 2 . Its f3 function 
reads 

-^=^ = T c, (3.13) 

giving 

c\t) = cl( 9 ^-)\ (3.14) 

Equation (|3.14|) shows that in the UV regime c 2 has a logarithmic RG running, governed 
by the coupling g. 

We expect that in any generic, weakly-coupled quantum field theory of Lifshitz type, 
including also theories with gauge fields and matter in D = 3, the running of c 2 will be 
qualitatively similar to (|3.14p . i.e. with a logarithmic dependence on the energy in the UV 
regime. 

3.2 Two scalar particles 

We will now show that theories with anisotropic scalings generically lead to different "speed 
of light" parameters (defined as coefficients of (dfa) 2 ) associated with different particles. 
More precisely, we will show that the RG evolution of the difference 5c 2 = c 2 — c 2 , even in 
the most optimistic case when Sc 2 = is an attractive fixed point, is generally too slow 
to give 5c 2 ~ with the needed accuracy. A severe fine-tuning in the UV for 6c 2 seems to 
be inevitable. 

In what follows we consider an extension of the single field model defined by the 
Lagrangian density (|3.ip to two fields fa and fa, imposing the Z 2 X Z 2 symmetry fa — > —(pi, 
i = 1,2. The Lagrangian is given by 

£ 2<t> = d + C 2 - gi2(fadifa){fadifa) - ^{d.fa) 2 ^ 2 - tj(di<h) 2 <fi ~ V 12 (fa,fa) , (3.15) 

where Ci ;2 are two copies of the Lagrangian appearing in (13. lj) for the fields fa and fa, 
and V± 2 (fa, fa) is an additional potential: 

Vu(fa,fa) = ^f<Pl4 + + k f^\ ■ (3.16) 

As can be seen, the Lagrangian contains a number of new interactions, which considerably 
complicate the analysis. In particular, the one-loop renormalization of the couplings k\, 
k 2 , k\ 2 and k 2 \ involves diagrams with all possible combinations of three insertions of the 



10 



quartic couplings gip, h\ %2 an d 512 as well as diagrams with one insertion of any of the 
order six terms and one insertion of any of the quartic couplings. Fortunately, at one-loop 
level, as in the single field model considered before, the renormalization of gip, h\ j2 , 512, 
c 2 and c 2 does not involve the k\, k 2 , &12, &21 couplings and therefore we do not need to 
compute the associated Feynman diagrams. In analogy to eq. (13.1ip . there will always be 
a choice of boundary conditions for the couplings at t = such that the model is stable 
all the way down to the far UV. 



Taking a 
we get 



a 2 for simplicity, after a straightforward but lengthy computation, 



092 



k 

8 
k 



(3gf + Ag 12 h 2 + hih 2 - 2h 



3g 2 + 4#i 2 /ti + h x h 2 



2h\ 



9i2 + hi(gi + 52) + hf + 2g 12 h, 



1,2 , 



0, 



912 



k 
16 



35 2 2 + 2^12(51 + g 2 ) + 3512(^1 +h 2 )-h 1 h 2 



(3.17) 



where all couplings have been rescaled by a factor 1/a 3 while keeping the same notation 
for the couplings for simplicity (i.e. we omit hats). The (5 functions for the c 2 couplings 
are easily computed to be 

Pe\ = ^ ( c i9i + 4hi) , 

Pel = l j(c 2 2 g 2 + clh 2 y (3.18) 

The RG equations (|3.17p do not admit a simple analytic solution in general. A class of 
exact solutions is however obtained by substituting the ansatz 

.9o 



gi(t) = g 2 (t) = g(t) -- 

hi(t) = h 2 (t) = h{t) 

/ 1\ 512,0 
512 it) 



1 — xkt ' 

hp 
1 — xkt 



1 — xkt 

in eqs. (|3.17p and solving the (now algebraic) equations for go>^o> 512,0 an d x: 

8xg - (3^o + 4^12,0^0 - hi) = , 
8xh - [gl 2 + 2g 12fi h + /i ( 2 5o + h)] = > 
16xgi2,o - [3^12,0 + 512,0(490 + 6/i ) - hi] = . 



(3.19) 



(3.20) 



In terms of go, ho and x, the running for c 2 



c\t) 



4 



go 



90+^0 

8x 



( c i + c i)/2 an d Sc 2 = (c 2 — c 2 )/2 is given by 

'g(ty 



5c\t) = 54 



go 



9Q~ h 

Sx 



(3.21) 
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The system ()3.20p is under-constrained (3 equations for 4 variables), so we fix one of the 
couplings, say go = 1, and look for solutions for x and for the other couplings /io,5i2,o- 
Taking go = r would simply rescale the solutions x — ► rx, ho — > rho, 512,0 — > r 9i2,o, so that 
the RG evolution of c 2 and 5c 2 is unaffected. A sufficient condition to get a (semi)positive 
definite interaction requires go, ho > and (512,0! < (50 + ho)/2. We find seven solutions, 
one of which is unstable and is disregarded. The remaining six solutions are 



1) x 

2) x 

3) x 

5) x 

6) x 



512,0 = . 



3 

4' 
5 

12' 

^(s^v 7 !?), 

^(77 + 5^17), 



512,0 



512,0 



512,0 



512,0 



= 1, 

1 

3' 



17), 
17), 



ho 
ho 
ho : 
ho 
h 



0; 

i; 
1 

3 ! 

:i(13=F3VTr); 

: 1(1 + VTT). 



(3.22) 



As can be seen from x, all six solutions correspond to couplings (and c 2 ) which grow in 
the UV in all cases. The deviation 5c 2 , instead, increases in the UV for 1), 3), 4) and 6), 
is constant in the case 2) and decreases in the UV in 5). Note that solution 1) reproduces 
the RG evolution (|3.1ip for a single field. 

We can also perturb the solutions found above and study the RG evolution of the 
fluctuations at the linear level. We focus on the case 1), since it is the only one giving rise 
to simple analytical results. We look at linear perturbations around the solutions putting 



5i = 5 + 5g + 8u , 
hi = h + 5h + 5v , 



52 = 5 - 5g + 5u , 
hi = h — 5h + 5v , 



and consider be 2 as a fluctuation. In this way, we get 
6g(t) Su(t) (g{t)' " 



<$5o 
6h(t) 
5h 

5c 2 (t) 



5uo 
5v(t) 



5o 
Sgn(t) 



5vo £512,0 

g(t) 



50 



5cl 



5o 



5o / 



512 



9M 
50 



512 + %2 , (3.23) 



1 



Sh 



• 



(3.24) 



From eqs. (|3.24p we see that the fixed-point is stable, with all fluctuations decreasing 
towards the IR. A similar study can be done for the other solutions, which are not all 
IR stable. From eq. (|3.24| ) we see that it is not enough to start with 5cq = at fj,o 3> 1 
to ensure that be 2 = near E ~ A. In addition, one has to ensure that also other 
perturbations around some fixed point are fine-tuned to zero. 
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(a) 



(b) 



Sc 2 



Sc 2 




t 



Figure 2: (a) RG evolution of Sc 2 as given by eq. (|3.24|) for eg = 1, Sg = 8h = 0, <5cq = 10 8 . (b) 
RG evolution of Sc 2 as given by eq. (|3.24| for c 2 — 1, 5g — Sh — — 1CT 2 , <5cq = 0. 

In section 5 we will argue that in the IR the RG evolution of 5c 2 , as given by the con- 
tributions of only weighted marginal couplings, essentially stops. Given the experimental 
bounds on 5c 2 for ordinary particles mentioned before, recovering Lorentz invariance in 
the IR with enough accuracy requires fine-tuning of 5c 2 , 5g and 5h in the UV to extremely 
small values. We illustrate this point in figures EJa) and[2^b), where the running of Sc 2 (t) 
over 40 orders of magnitude for an arbitrary given choice of boundary conditions is shown. 

3.3 Case of particle with no self-interactions 

The two-scalar model can be adapted to the case where one of the particles, say (j>\, has no 
self-interaction terms, i.e. g\ = 0. This case represents a situation of physical interest: in 
electrodynamics the photon has no self-interaction term and one may wonder if the speed 
of light will still significantly depend on the energy. Here we will show that the energy 
dependence of either speed-of-light parameters c\ or C2 does not rely on the presence or 
absence of self-interaction terms. 

Consider the beta functions (|3.17j) . An exact solution can be found by setting 



gi = h 2 = gn = . 



(3.25) 



We are left with two equations 






(3.26) 



Next, we write p = fh\, so that the second equation takes the form 





(3.27) 
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One obvious solution is / = which requires / = /+ or / = /_ . The solution with / = /_ 
is not physical, since / = 1 + gijh\ > 1- Setting / = /+ leads to a solution similar to the 
class of solutions considered in (|3.19l) . 

,.\ 9o , M h f + h 

92{t) = T^Ut> hl ® = T^kt> x = — - (3 ' 28) 

Now consider solutions where / is not constant. Writing 

1 = ^ = 4--^!, (3.29) 
J dhi dhi 8 1J y ' 



(3.30) 



eq. (|3.27p becomes 

f df _ 2 

(/-/+)(/-/-) dfei fci ' 
We can now integrate this equation and obtain 

!+ 

( f _ f \2(/+-/_) 

/»i = - J + > f , (3.31) 

(/ _ / _ )2 (/ + -/_) 

where &o is an integration constant. Inserting eq. (|3.31[ ) into the first equation in (|3.26|) 
gives a first order differential equation for / that can be integrated. The result is 

kk f°° df 

Jf (/'-/+) + 2 (/ + -/-)(/'_/_ ) 1 m=J=) 

We have chosen as the ultra- high energy scale at which / — > oo, so that / is defined in 
the interval — oo < t < 0. The integral (|3,32p defines / = f(t) and hence h\{t) = hx[f(t)] 
and g2(t) = (f(t) — l)hi(t). It can be expressed in terms of hypergeometric functions, 
but its explicit expression is not needed in order to see the qualitative behavior of the 
couplings in the regime i C 0. In this limit, / — ► f+ and hence one approaches the 
constant / solution (|3.28|) discussed above. For large (— t), one has the behavior, 

(f - f + )W&=l „ hl „ to „ * (3.33) 



Hence both h\ and 52 decrease towards the IR. 

Next, we consider the RG equations (|3.18p for c\ and C2, which simplify to 

(3.34) 

Integrating these equations using the behavior of 52 and hi given in eq. (|3.33p . one can 
see that C2 goes to zero like l/|i| fc , with k > 0, while cf approaches a constant value for 
t < 0. 
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4 UV behavior: a UV free model with z = 2, D = 10 

In this section we look for a weighted renormalizable scalar field theory which is UV 
free. In order to simplify our analysis, we consider a (j) 3 model which is weighted-counting 
renormalizable in D = 10 spatial dimensions. Although the (f) 3 model is manifestly an 
unphysical theory, having an unbounded potential, this instability does not affect the RG 
equations, which can then be formally studied and hopefully seen as a toy laboratory 
for more complicated stable UV free theories, such as Yang-Mills (YM) theories. As in 
section 3, we first study the one-loop RG evolution of a single field, which will enable us 
to compute the RG evolution of c 2 . Then we shall consider a two-field model, where 8c 2 
will also be considered. 

4.1 One scalar 

The most general renormalizable Lagrangian reads 



As in section 2, we use a dimensional regularization in the spatial directions only. We will 
add counterterms to subtract the poles in 1/e only, with no finite term. We study the 
theory in the energy range E 3> 1, where the analysis is reliable. 

The wave-function renormalization in the <f> 3 theory is non-trivial and hence the rel- 
evant one-loop graphs to compute are those associated with the two and three point 
functions. These graphs are computed exactly along the same lines followed in section 3, 
so we will just report the results for the (3 functions and the anomalous dimension 7 of (p. 



(4.1) 



We find 7 = /i A 2 /(64o 5 ) and 

A 2 A 2 
(3 a 2 = 2u a —a 2 , p x = -^A^-A , 



(4.2) 




(4.3) 



X = —uj x x 



2 



(4.4) 



with u). 



•x — 



2u)x + 5u> a , and thus 



x(t) 



1 + X(jU x t 



(4.5) 



Plugging eq. (|4.5p into eqs. (|4.2p give 




(4.6) 
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The coupling a 2 increases in the UV, while A and the effective coupling x are UV free. 
Let us now study the evolution of c 2 . Its associated (3 is 



where u> c = 3Zio/16. Hence 



P c 2 = U) C XC , 



(4.7) 



(4.8) 



The RG evolution of c 2 is again governed by the weighted marginal couplings of the theory. 
Interestingly enough, c 2 logarithmically increases in the UV, as in the model of section 3. 



4.2 Two scalars 

We can add a further scalar rj to the <f> 3 model and, for simplicity, we impose a discrete 
Z 2 symmetry under which rj — > — rj. The Lagrangian is given by 

I ^2 „2 ^2 



C = 



^ 2 -f(Ar?) 2 -f(^ 



7U2 



rn 



-rj 



A .3 

" 3!^ + 
-rj (p. 



(4.9) 



2 ' 2 

The computation of the (5 functions for A, A, a 2 , a 2 , c 2 and c 2 is straightforward, although 
a bit laborious, mainly because we keep a 2 / a 2 in general. We find 

1 



fa 
fa 



-X 3 F(a, a) + A 3 ( -F(a, a) + ^K(a, a)) + ^XX 2 K(a, a) 



-XX 2 F(a, a) + -X XK(a, a) + A d -F(a, a) + 2K(a, a) + -K(a, a 



A 2 (a 2 if(a, a) + Q(a, a)) + A 2 (a 2 K(a, a) + Q(a, a)) 

2 ' 



4 ' 



= A 2 (a 2 K(a,a) + Q(a,a)) 



fa = 
fa = 



A 2 (c 2 K(a, a) + R(a, a, c, c)) + A 2 (c 2 iT(a, a) + i?(a, a, c, c)) 
A 2 (c 2 K (a, a) + R(a, a, c, c)) — , 



no 



where 



F(a 1 ,a 2 ) 
K(a 1 ,a 2 ) 
Q(a 1 ,a 2 ) 
R(a 1 ,a 2 ,c 1 ,c 2 ) 



2ai + a 2 
a\a 2 (ai + a 2 ) 2 ' 
1 

a\a 2 {ai + a 2 ) 3 ' 

(af + 5af a 2 + lOaf a 2 , + Said 3 , + a|) 

5aia 2 (ai + a 2 ) 5 
cfa^a? + 12afa 2 + 8aia| + 2a|) + (1 <-> 2) 
5afa^(ai + a 2 ) 4 



^10 
T 



(4.10) 



(4.11) 
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When A = 0, eqs. (|4.10p collapse to eqs. (|4.2p and (|4.7p . as expected. Finding the general, 
analytic solutions to eqs. (|4.10p is complicated. Interestingly, eqs. (|4.10p admit, neverthe- 
less, a simple fixed-point solution given by A = A = A(i), a 2 = a 2 = a 2 (t), c 2 = c 2 = c 2 (t). 
The RG equations for A, a 2 and c 2 are precisely given by eqs. (]4.2p and (|4.7j) . provided one 
makes the substitution u\ — > 2u;a, w a — ► 2u; a , and u; c — > 2a> c . For completeness, below we 
report the corresponding solutions: 

Mt)=um w , ? {t)= 4mY w . (4.i2) 

where 

x(t) ^^1= p ■ (4.13) 

The fixed-point solution (|4.12p is unstable under small deformations. In order to show 
that, we consider the following "symmetric" perturbations 

A = A + 5A, A = A-5A, 

2 -2 , i 2 ~2 -2 s 2 

a = a + oa , a = a — da , 

c 2 = c 2 + 5c 2 , c 2 = c 2 -5c 2 , (4.14) 

keeping up to linear terms in the perturbations. We get 

A 2 

Psx = -2w 5A — SX, 
a 

P8ai = ^^ma 2 5\ + W5a 2 \), 
loO a 

(3 5c 2 = l ^(6c 2 5X + 5c 2 X), (4.15) 
16 a 

where ujg\ = 7iio/64. The UV evolution of the couplings is therefore given by 



SX(t) = 5X 



x 



alSX \ 


(m 


Ao / 




(%6X \ 


f x (ty 


Ao / 





4 2r 62 



Ao V x o 

4 28 

5c\t) = (5c 2 -^\r-^) ~^ + <Bo(my^ (416) 

V A / V / A V / 

The fixed-point is stable in the UV only along SX. The situation is similar to the one 
encountered in the model of section 3. The RG evolution of 5c 2 in the UV does not help 
in alleviating the fine-tuning needed to get Sc 2 small enough for energies near A. 

We expect that the general lessons for the behavior of the different couplings under the 
Renormalization Group learned in this section should be qualitatively shared by a large 
class of UV free quantum field theories, including the relevant case of perturbative YM 
theories coupled to matter. 
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5 IR behavior of Lifshitz-type theories 



As well known, in conventional Lorentz invariant theories, the MS scheme has to be used 
with care in studying the evolution of the couplings in presence of mass terms when 
E <C to, since the decoupling of massive particles is not manifest. Indeed, the MS (5- 
functions, being mass-independent, are formally the same for any E, while for E <C to the 
physical coupling does not effectively run anymore. In conventional theories, a reasonable 
approximation to overcome this problem and still use the MS scheme is to take an effective 
approach, where for E <C to the heavy particle is integrated out and its contribution to 
the running neglected. One could also make a more refined study of the transient region 
around E ~ to (we will not do it in our similar case below). Since we are investigating 
here unconventional theories of the Lifshitz type, it is instead preferable to use a more 
physical renormalization scheme, such as the momentum subtraction scheme, where the 
/^-functions are mass dependent and the decoupling is manifest, even if the associated 
one-loop computations become necessarily more involved. In the present case, because 
the propagator is given by 

^-c^-^-to 2 )" 1 , (5.1) 

IR effects are expected as soon as c 2 A 2 k 2 ~ a 2 k 4 , i.e. at momentum scales of order cA/a. 
The scale cA/a plays now the role of to. As explained below, the underlying reason is 
that the term c 2 k 2 in the propagator itself provides an IR regularization of the amplitudes 
even when to = 0. 

5.1 (3- function in the IR in conventional (ft 4 theory 

Let us begin by briefly reviewing the standard computation of the one-loop /^-function of 
the Lorentz-invariant theory in 3 + 1 dimensions in the momentum subtraction scheme. 
Let I^ 4 ) be the tree+one-loop+counterterms 1PI four-point function. The physical cou- 
pling A can be defined as the value of at some given energy scale. In terms of the 
Mandelstam variables s, t, u, we can define 

r< 4 > ( s = t = u = -^ 2 ) = (n) = -x. (5.2) 

Here we used the standard trick of introducing the renormalization group scale at some eu- 
clidean value of the kinematical invariants to circumvent threshold singularities. Eq. (|5.2p 
fixes the value of the counterterms, so that 

T^(s,t,u) = -A + r[ 4) ( S ,t,n) - if V), (5.3) 

where is purely the 1-loop contribution to . Since (3 m i and the anomalous dimen- 
sion of 4> vanish at one-loop order, being T) momentum-independent, inserting eq. (|5.3j) 
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in the Callan-Symanzik (CS) equation satisfied by I^ 4 ) gives directly (3\ = fidX/dfj,: 



8r w 



Using standard techniques, we get 



; - fi >2 / ^kE_ [\ /i 2 x(l-x) [\ y 2 x(l-x) 

A ' (2vr) 4 7o + m 2 + ^ 2 x(l - x)] 3 16vr 2 y o m 2 + /x 2 x(l - x) ' 1 j 



The integral in x in eq. (|5.5j) can analytically be performed, but it is not necessary to 
do so in order to see the UV and IR behavior of f3 X - For /i 2 > m 2 , f3 x ^ 3A 2 /(16tt 2 ), 
in agreement with what one would have obtained with, say, a computation in the MS 
scheme. For /i 2 <C m 2 , f3\ ~ A 2 /(327r 2 )(// 2 /m 2 ). The behavior of A below m, its freezing, 
can essentially be understood by noticing that m acts as an IR regulator to the one-loop 
graph, which would otherwise have an IR divergence when /i — ► 0. 

5.2 Computation of /3-functions of Lifshitz-type theories in the IR 

As we have already mentioned at the beginning of section 5, even in the absence of mass 
terms, the RG evolution of coupling constants presents at least two regimes, characterized 
by fi S> c/a (UV) and fi <C c/a (IR)EI Here we consider the (ft 3 model of section 4 in detail 
as an illustrative case. As it will shortly be clear, the important qualitative aspects of the 
results should be completely general and apply to any Lifshitz-like theory, including the 
model in section 3. We define the renormalized field <j) an d the couplings A, a 2 and c 2 as 
follows: 



ar( 2 ) 



d(p' 



,0\2 



_ ] i d 4 r< 2 ) 

~ ' 4! dp 4 



-a 2 



i a 2 r( 2 ) 







2 dp 



,2 



-c 2 



r( 3 )[(rf) 2 = -a;(/i 2 ),P2,3 =Pi,2,3 =0] = r( 3 V) = -A, (5.6) 

where r( 2 ) and T^ 3 ) are the tree+one-loop+counterterms 1PI two- and three-point func- 
tions and stands for the subtraction point (p ) 2 = —u>(fi 2 ) = — a 2 /i 4 — c 2 /i 2 , p = \p\ = 0. 
We define the subtraction point at vanishing spatial momenta for technical reasons, since 
the computation is greatly simplified in this way. Due to the Lifshitz nature of the theory, 
the energy p° should depend quadratically (for z = 2) on the sliding scale /i, so in the UV 
we get uj(fj?) ~ a 2 /x 4 . The factor a 2 has been inserted just to slightly simplify the analysis 
that will follow, but it is by no means necessary. In the IR we get to(/J- 2 ) — c 2 ^ 2 . In 
both cases, we have chosen the subtraction point at an euclidean energy scale to circum- 
vent threshold singularities; this is a non-relativistic analog of the more familiar condition 
(|5.2p . Using the CS equations satisfied by I^ 2 ) and we can derive a 2, [5 C 2 and 



:> In units of A. In the following, A will be restored in some formulas when convenient. 
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the anomalous dimension 7 of eft in terms of the purely one-loop contributions (denoted 
by rj 2) and rj 3) ) to and . We have 

1 <cJ R OT S' 

1 dr (2) "" 1 dr (2) " 

^2 = a + 27a , /3 C 2 = ~2^—^- + 2 ^ c • ( 5 - 7 ) 

To simplify the notation, here we have denoted by a dot and a prime a derivative with 
respect to (jP) 2 and p, respectively. We now show the computation of 7 in some detail. 
One has 

1 ,,2,.J(, ,2\J2(-\ _ „.\2 



d W kdk E f\ fl 2 Lo'(^ 2 ) X 2 (l 



7 = 3A / / dx 



■r 



(2vr) n J [k% + a 2 k 4 + c 2 k 2 + cu(fi 2 )x(l - x)f 

15 ''° A2 '^r^^M^^r. ,5.8) 



32 J J [a 2 k 4 + c 2 k 2 + w(p 2 )x(l - x)} 7 / 2 • 
Computing the integral in x, we find 



3C- 



7 = A 2 / 10 / dk - fcVW) . (5.9) 

Note that the integrand in eq. (|5.9p is UV and IR finite for any value of fi 2 > 0. In the UV, 
fi S> cA/a, the c 2 terms can be neglected, in which case the integral is easily performed 
giving 

7(UV) "1^' » y> -a A - (5 ' 10) 
In the IR, one can neglect the fi 2 term appearing in the denominator of the integrand of 
eq. (|5.9p . in which case we get 

Equation (|5.1ip implies that the RG evolution of 7 in the IR is essentially frozen, in 
complete analogy to what happens in standard Lorentz invariant theories below the mass 
scale. Equation (|5.1ip is best understood by noticing that eq . (|5.8j) is well-behaved in the 
IR, due to the c 2 k 2 term that acts as an IR regulator and forbids the presence of any IR 
singularity for k — » 0. Like in the usual Lorentz-invariant </> 4 theory considered before, 
the IR finiteness is responsible for the freezing of the coupling. The essential aspects 
of the above results should be general and apply to any weighted marginal operator in 
Lifshitz-type theories. 

A computation very similar to the one above (but algebraically more involved) gives 



3 m~ k , A " »\ 5 (iR)^ 2fc g AV 2 (IR)„ 2fc c A 2 ji 2 
Vc 2 A 2 ' Pa2 ~ cfic 2 A 2 ' Pc2 ~ a 3 c 2 A 2 



2 
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where k x = -7l 10 /480, k a = 183/i /22400 and k c = -3Zi /400. We have checked that 
/3^ UV \ /3^2 JV ' > and /3 2^ completely agree with those found in the MS scheme in section 
4, eqs, (!4.2p and (|4.7j) . Notice that it is crucial to take p° oc /i 2 in the UV regime to get 
agreement with the MS scheme. By taking p° oc fi, we would have obtained a mismatch 
by a factor 1/2 in the /3^ uv ' -functions. The effective IR coupling is 

A 2 

V = > ( 5 - 13 ) 

arc? 

to be compared with the UV effective coupling x = A 2 /a 5 . Contrary to x, however, y is 
not well-defined, in the sense that it sensitively depends on the choice of subtraction point, 
i.e. the choice of [i. This is obvious, considering that in the IR the /3-functions present 
an explicit dependence on fi. It is straightforward to solve the RG equations in the IR. 
Restoring A, we get 



y(A) 



X(jj) = A(A) 



a 2 {n) = a 2 (A) 



l + ^y(A)(l-^)/2 

v(p)' 



y(A) 
y(A) 



c 2 (,) = c 2 (A)(M)\ (5.14) 

where k y = 2k\ — 2k c — 3k a . Eqs. (|5.14p exhibit the fact that in the IR regime the RG 
evolution of a 2 , c 2 and A, as induced by A itself, freezes. More precisely, for fi <C cA/a - 
modulo a small threshold effect shifting the IR values of the couplings from the UV ones 
- the energy-dependence of the low-energy couplings is proportional to (/u/A) 2 . 

5.3 IR effects of relevant couplings 

In the previous subsection we have shown that weighted marginal operators become inop- 
erative for n <C cA/a. This is in some sense expected, since in the IR one should recover 
the usual classification of operators in terms of canonical rather than weighted dimensions. 
What is weighted marginal in the UV becomes standard irrelevant in the IR. However, 
care has to be paid to the weighted relevant operators which become standard marginal 
in the IR. These operators deserve a separate discussion. They are obviously negligible 
in the UV but they can efficiently transmit the UV Lorentz violation to the IR theory. 
Indeed, as we will see, Lorentz-symmetry-breaking weighted relevant operators which be- 
come standard marginal in the IR will in general mediate Lorentz violation effects to the 
whole IR Lagrangian. 
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In order to illustrate the effect, we consider a simple IR toy model in 3+1 dimensions, 
consisting of a fermion interacting with a scalar by means of a Yukawa interaction. We 
imagine that, say, the scalar is described at high energy by a Lifshitz-like dynamics with 
some z > 1. The Yukawa coupling is standard marginal in the usual sense, but weighted 
relevant. We want to study its effect in the IR where, as we have just shown, the weighted 
marginal couplings of the Lifshitz theory have a negligible beta function. We assume that 
the only remnant of the non-Lorentz invariance of the UV completed theory is a tiny 
difference in the speed of light of the fermion and scalar: cw, = 1, = 1 + 5c, with 5c <C 1. 
The Lagrangian is simply 

1 c 2 

£ = 2^) 2 " y W) 2 + - hH) ~ (5-15) 

We will now show how, due to a non-vanishing value of 5c, a logarithmic running in 
cv, is induced. We work at linear order in 5c and in the MS scheme. We omit details, 
since the one-loop graphs of this model are straightforward and can be found in standard 
textbooks. We denote by and Z^ the wave-function renormalization constants of <p 
and tp, associated with the renormalization of ((ft) 2 and ^00^, by Z C( p and Z c ^ the wave- 
function renormalization constants associated with (di4>) 2 and -ip0ip and by Z g the vertex 
renormalization constant. We get 

-2 



Z c ^ = l- J l-^\( 2 -5c + 0(5c 2 : 



Z c ^ = l-j£^±(85c + 0(5S) 

~2 



9_ 

(4tt) 2 e 



Z g = l + 7 f—-[2 + 0{5c)). (5.16) 



From eqs. (|5.16p we get 



0^ ^9^ ^9^ 

0c, = -^Sc, Psc = ^8c, (3 g = ^- (5.17) 

Defining a = g 2 /{Att), we finally get the following RG evolution for c^(t), 5c(t) and a(t): 

1 - ^a t 



5c(t) = 5cq 



a(t) 
a 



C^[t) — C0 iO — 



a(t) 



7_ 

15 

- 1 



(5.18) 



a 

Equations (|5.18|) shows that any small Lorentz symmetry breaking term in the IR theory 
(coming from the underlying UV theory) induces, by quantum effects, an energy-dependent 
speed of light for all particles sensible to the breaking term. From this example it should 
also be clear that the effect is general. 
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6 Conclusions 



In this work we have carried out the analysis of the one-loop RG evolution of Lifshitz-like 
theories, by mainly focusing on two specific Lorentz violating scalar field theories. Our 
primary interest was to look for possible simple mechanisms (see also |15| ) which would 
alleviate the fine-tuning otherwise needed in these theories to recover Lorentz symmetry 
at low energies with the needed accuracy. We have focused our attention on two particular 
operators, the spatial kinetic terms c 2 2 (<90i j 2) 2 . There are essentially two regimes of in- 
terest. In the UV Lifshitz-regime, c 2 and 5c 2 logarithmically run with the energy, and the 
running is governed by weighted marginal operators. In the IR, if all sources of Lorentz 
symmetry breaking are vanishing, one recovers Lorentz invariance, which clearly forbids 
any scale-dependence for c 2 . However, due to the effects of standard marginal couplings, 
any small Lorentz symmetry breaking term in the theory leads to phenomenologically 
unacceptable logarithmic dependence of c 2 on the energy scale. As it has already been 
pointed out, our considerations seem to be very general and independent of additional 
structure and/or symmetry that would be present in more realistic theories, such as gauge 
symmetries. It should be clear that, although we have been focusing on the spatial ki- 
netic term operator, our considerations may be equally applied to other Lorentz violating 
parameters associated with other operators in a theory (this could lead to the need of 
additional fine-tuning). 

In section 5.3 we computed the evolution of c in the IR in a simple low-energy scalar- 
fermion theory, assuming the presence of a small 5c, coming from the non-Lorentz invariant 
UV theory. In this example 5c is driven to small values at low energies, see eq. (|5.18p . 
A logical possibility is that in the Standard Model the RG flow could similarly drive the 
5cij = Ci — Cj corresponding to any i, j pair of particles to sufficiently small values, below 
the existing phenomenological bounds. Although not completely excluded by our analysis 
(since we are not computing RG in a Lifshitz-type Standard Model), this possibility seems 
unlikely, since the logarithmic running of 5c at low energies is too slow to drive a 5c of order 
1 at the scale A to values compatible with the stringent existing bounds 5c < 10~( 21 ~ 23 ). 
Moreover, in a UV free theory, the same mechanism of section 5.3 would give a 5c that 
increases towards the IR. 

An order of magnitude of other potential phenomenological problems can also be ob- 
tained by applying our considerations to the photon in a weighted renormalizable version 
of QED. At low energies we can neglect higher orders terms in p (always present in Lifshitz 
type theories), so that the photon dispersion relation will be typically of the form 



,2 



c"(t)p 2 




where oj is the energy and 




(6.2) 
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with some unspecified constants / and r. The behavior eq. (|6.2p should be induced by one 
loop effects as in section 5.3 once a small 5c between the photon and a charged particle is 
introduced in the Largrangian. This behavior can be compared with known experimental 
constraints. An interesting constraint on the energy dependence of c 2 at moderately high 
energies has been recently given by the FERMI experiment, detecting the photon spectrum 
of the gamma ray burst GRB 080916C at red-shift z = 4.35 [130. By measuring the time 
delay (~ 10 sec) between the "low-energy" (< 1 MeV) component of the burst with respect 
to the high energy one (~ 10 GeV), we get the following rough experimental boundgj 

|c 2 (lMeV) - c 2 (10GeV)| < 10~ 17 . (6.3) 

By taking r ~ 1 in eq. (|6.2j) . eq. (|6.3p approximately gives the constraint on / described in 
eq. (jl.4p . Similar (although milder) bounds exist for other particles as well. 

One might think that, in the context of the Hofava theory of gravity, the bounds we are 
finding can be evaded for all the Standard Model (SM) particles by simply not introducing 
the higher derivative spatial interactions responsible for the effects studied in this paper. 
Being the SM renormalizable, there is no real need of introducing them. Aside from the 
lack of a clear principle besides this choice, radiative effects induced by graviton loops 
should nevertheless generate these Lorentz-violating couplings in the SM sector at some 
loop order. 

Finally, let us conclude by noticing that the one-loop corrections to c 2 we considered 
are quadratically divergent. If new physics above the scale A is assumed, this can give rise 
to a naturalness problem, similar to the standard one affecting the Higgs boson mass. 
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